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Abstract 

The purpose of this paper is to investigate the limiting distribution functions for a polynuclear growth 
model with two external sources, which was considered by Prahofer and Spohn in |l3| . Depending on the 
strength of the sources, the limiting distribution functions are either the Tracy- Widom functions of random 
matrix theory, or a new explicit function which has the special property that its mean is zero. Moreover, 
we obtain transition functions between pairs of the above distribution functions in suitably scaled limits. 
There are also similar results for a discrete totally asymmetric exclusion process. 
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1 Introduction 

Our main object of study is the following combinatorial problem. Fix three real parameters t > and a± > 0. 
We construct a random set of points in the unit square [0, 1] x [0, 1] C M 2 , as follows. Let P(A) denote a Poisson 
variable of density A. We select P(t 2 ) points at random inside the square (0, 1) x (0, 1), P(a + t) points at random 
on the open bottom edge (0, 1) x {0}, and P{a~t) points at random on the open left edge {0} x (0, 1). Hence 
no point is selected from the lower left vertex and from the closed top and right edges. For example, in Figure 
[j], 5 points are selected inside the square, and 2 and 1 points in the bottom and the left edges, respectively. A 
(weakly) up/right path is given by a sequence of points such that each point is (weakly) above and to the right 
of its predecessor; thus the solid line in the figure is a (weakly) up/right path from the lower left vertex to the 
upper right vertex. The length of a (weakly) up/right path is defined by the number of points in the path; thus 
the up/right path of the example has length 4. 

We define L(t) be the length of the longest (weakly) up/right path in this random configuration of points. 
(This, of course, implicitly depends on the values of a±.) E.g., the solid line in Figure |l| is the longest up/right 
path of the example. In general, there may be more than one longest path, but we are only interested in the 
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Figure 1: A points selection and the longest up/right path 



length L(t), hence the degeneracy is not an issue. We note that the above process can be thought of in an 
alternative way as a Poisson process of intensity 1 in the open , together with a Poisson process of intensity 
a+ on the open half- line R+ x {0} and a Poisson process of intensity a_ on the open half-line {0} x R_. Then 
L(t) is equal to the length of the longest (weakly) up/right path from (0,0) to (t,t). 

The main interest in this paper is the statistics of L(t) as t — > oo, as this problem arose in a polynuclear 
growth (PNG) model considered by Prahofer and Spohn in PNG is a simplified model for layer by layer 
growth in one spatial dimension. At each random nucleation position, an island of height 1 is formed and 
spreads laterally with speed 1. When two islands meet, they form one island and keep spreading with the same 
speed. The question is the fluctuations of the height in the large time limit. 

As a special case, suppose that a single island starts spreading at the origin and that further nucleation 
takes place only on top of this ground layer. Furthermore, suppose that there are external nucleation sources 
at the two ends of the ground layer. Then the height, for example, at the origin at time y/2t is equal to L(t) 
in the above point selection process. We refer the interested readers to the papers |l3|, [l2| for details of the 
mapping from PNG to L(t) and many other related works. 

In the above point selection process, a special case is when a± = : P(t 2 ) points are selected inside the 
square with no points on the edges. (This corresponds to a PNG model without external nucleation sources.) 
It is interesting that in this case, there is a combinatorial interpretation. Consider a random permutation of 
Sjv- If we take N as Poisson variable P(t 2 ), the longest increasing subsequence of a random permutation of Sn 
has the same statistics as L(t) (see Q and references therein). The limiting fluctuation of L(t) in this case is 
obtained in g] : there is a distribution function Fque^) such that 

limp|^ < x) = F GUE (x), a± = 0. (1.1) 

The convergences of all the moments are also proved in the same paper. Here Fque is the so-called GUE 
Tracy- Widom distribution function of random matrix theory. See Section || for the explicit formula of Fque, 



equation (2.9) 



The "GUE" in Fque refers to the Gaussian unitary ensemble, the set of all N x N Hermitian matrices 



2 



together with the probability measure defined by 



— e- trM2 dM 



(1.2) 



where Zm is a normalization constant (see e.g. flfj| , |j|). Tracy and Widom in 



|jl4| proved that as N — > oo, the 



properly centered and scaled largest eigenvalue of a random GUE matrix converges in distribution to Fgue. 



random GUE matrix have the same fluctuations. 

In this paper, we obtain similar limiting distributions in the presence of points on the edges. When a± > 0, 
the longest path begins by following one of the edges for a while, then enters the square. Hence when a± are 
small, the effect coming from edges is small, and we expect L(t) to have the same statistics as when a± = 0, 
i.e. GUE fluctuation. But when ot± is large, the longest path lies mostly on one of the edges, hence we expect 
Gaussian fluctuation. It turns out that the critical case is when a± = 1. We need to distinguish four different 
cases to state the results. In each case, a± are fixed. 

• When a± < 1, we obtain Fgue m the limit. 

• When either of a± is greater than 1, we obtain Gaussian fluctuation. 

• When one of a± is equal to 1 and the other is strictly less than 1, we obtain Fqoe(^) 2 - In the above 
definition of the GUE, if we replace Hermitian matrices by real symmetric matrices, we obtain the Gaus- 
sian orthogonal ensemble (GOE). The limiting distribution of the (properly centered and scaled) largest 
eigenvalue of a random GOE matrix is given by Fqoe in |1| (see Section || for the explicit formula). The 
above limit Fqoe(^) 2 can be interpreted as follows : take two random GOE matrices and superimpose 
their eigenvalues. The largest of the superimposition of eigenvalues has limiting fluctuation Fgoe 2 - 

• When a± — 1, we have a new limiting distribution which we denote F , for which we do not yet have a 
random matrix interpretation. See Section || for the definition and discussions. 

More explicitly, we have the following theorem. 

Theorem 1.1. For each fixed a±, as t — ► oo, we have the following results, 
(i). When < a± < 1, 



Therefore the above result (1.1) implies that in the limit, L(t) with a± = and the largest eigenvalue of a 




Fgue (a;), 
FgoeM 2 , 



< a ± < 1, 



a± = I. 



aq_ = 1, < a_ < 1, or a_ = 1, < aq_ < 1, 



(1.3) 



(ii). When at least one of a± is greater than 1, setting a — max{a + ,a_}, 




(1.4) 
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Remark. As a tends to oo, L{t) tends to a variable with mean and variance at] i.e., to P(at), corresponding to 
the number of points on the interval (0, 1). For finite a, the longest path leaves the interval near 1 — a~ 2 , then 
includes roughly 2a~ 1 t points inside the square, thus giving L(t) ~ P((a — of agreeing with Q , 

The functions, FquEj Fqqe and F , in the above theorem are defined in Section |^, and we discuss their 
properties in the same section. 

The above theorem shows that there are certain transitions around the points a± = 1. It is of interest to 
investigate these transitions in detail. In Section |[ under proper scalings of a± — ► 1, we obtain new classes of 
distribution functions interpolating the functions of the above theorem. 

In addition to the above point selection process, there is a closely related lattice directed polymer model in 
2-dimensional space, which is a generalization of the model considered by Johansson in j^]. This model and a 
related exclusion process are discussed in Section ||. 

Finally the proofs of the theorems are discussed in Section [| 

Notations. In many papers, the functions Fque and Fqqe are denoted by F 2 and F%, respectively. 

Acknowledgments. The authors greatly appreciate many enlightening conversations and communications 
with Michael Prahofer and Herbert Spohn who brought our interest to this problem. 

2 Limiting distribution functions 

Let u(x) be the solution to the Painleve II (PII) equation, 

u xx = 2u 3 + xu, (2.1) 

with the boundary condition 

u(x) ~ —Ai(x) as x — + +oo, (2-2) 

where Ai is the Airy function. The proofs of the (global) existence and the uniqueness of the solution were 
established in ||. The asymptotics as x — > — oo are given by (see e.g. []|, fjj) 

- e -(4/3)- 3/2 

u{x) = -\j-^{^ + °{^2^ asi->-oo. (2.4) 

— (2/3)^^/^ 

Recall that Ai(x) <~ e 2 ^ x i/4 asu +oo. Define 



u(x) — — Ai(x)+0[ YJZ )' as a; — > +oo, (2.3) 



v(x) := / {u{s)) 2 ds, (2.5) 

J oo 

so that v'(x) = (u(x)) 2 . We note another expression 

v(x) = u{xf + xu{xf - {u'{x)) 2 (2.6) 
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which can be obtained by noting that the difference (i) has derivatives equal to zero by the Pll equation, and 
(ii) becomes zero as x — > oo by ( |2.3| ). 

The Tracy- Widom distribution functions are defined in terms of u and v. 



Definition 1 (TW distribution functions). Set 

1 



F(x) 
E{x) 



exp 
exp 



v(s)ds 

i 

u(s)ds 



exp 



(s — x)(u(s)) 2 ds 



and set 



Fgue(z) := 
Fgoe(z) := 



F(x) 2 
F(x)E(x) 



exp 



(s — x)(u(s)) ds I , 



(F 2 (x)) 1/2 e^™ u ^ ds 



GSE 



(x) 



F(x){E(x)- 1 + E{x)}/2 = (F 2 (x)) 



1/2 



/2- 



(2.7) 
(2.8) 

(2.9) 

(2.10) 
(2.11) 



It is proved by Tracy and Widom in |l4|, [l5| that under proper centering and scaling, the distribution of 
the largest eigenvalue of a random GUE/GOE/GSE matrix converges to Fgue(^) / Fqoe(^) / Fgse(£) as the 
size of the matrix becomes large. The readers are referred to [ [To| |B| for definitions of various random matrix 
ensembles and their basic properties. We note that from the asymptotics (2.3) and (2.4), for some positive 
constant c, 



F(x) = l + 0(e~ c: 

E{x) = l + 0(e- c: 

F(x) = 0(e- c ^ 3 ) 

E(x) = 0(e- c W 3/ '' 



.3/2, 



as x — > +oo, 
as x — > +oo, 
as x — * — oo, 
as x — > — oo. 



(2.12) 
(2.13) 
(2.14) 
(2.15) 



Hence in particular, as x — > +00, all the above three functions become 1, and as x — > —00, they become 0. 
Monotonicity follows from the fact that they are limits of sequences of distribution functions, and therefore 
(2.9)-( 2.11 ) are indeed distribution functions. 



We need a new distribution function for the case when a± — 1 in Theorem 1.1 
Definition 2. Set 



F (x) = {l-(x + 2u\x) + 2u(x) 2 )v{x)}(E{x)) 4 F GVE {x). 



(2.16) 



The asymptotics (2.3)-(2.5) imply that Fo(x) has limit 1 as x — > +00 and as x — > —00. The monotonicity 
of Fo follows from the fact that it is a limit of distribution functions in Theorem [O]. It would of interest to have 
random matrix interpretation of the function Fo, but so far we have been unable to identity Fo as a quantity 
arising in random matrix theory. 

One special property of this distribution function is that it has mean zero. 
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Proposition 2.1. We have 

/oo 

xdF (x)=Q. (2.17) 
-oo 

Proof. We note that the term in front of v(x) in the definition of F has another expression : 

x + 2m' (z) + 2u(x) 2 = E(x)- 4 ( E(t) 4 dt. (2.18) 

This follows by noting that y(x) := x + 2u' 4- 2w 2 satisfies 

2/'(x) = l + 2u(x)y(x), y(x) = -^=(l+o(l)), x -> -oo. (2.19) 

V— 2a; 



Then we have 



F " ( -' 1 = H F gve(x)J E(t) A dt^, (2.20) 



which upon integrating gives 

F (t)dt = F GUE {x) I E(t) i dt = F GVE (x)E(x) 4 (x + 2u'{x)+2u{x) 2 ). (2.21) 



Since ird Fo (ir) = lim x _>oo [x Fo(ic) — Fo(y)dyj from integration by parts, subtracting ( 2.2l| ) from 
xFq(x) and taking the limit x — > oo, we find that Fo has mean 0, as required. □ 

Remark. The mean zero property of Fo(x) was suggested in Jl^| by numerical computation. Moreover, by an 
indirect argument for PNG with a± = 1, it is shown that the average of L(t) is exactly 2t (see 11 ), which 
implies that Fo has mean zero. We note that the means of Fque and Fqoe are —1.77109 • • • and —0.76007 • • • , 
respectively. 



3 Around the transition : a± — > 1 

In this section, we investigate the transition around a± = 1 in detail. 

To state the results, Theorem 3.3, we need some preliminary definitions. Let F be the real line R, oriented 
from +00 to —00. Let m(-;x) be the solution to the Painleve II Riemann-Hilbert problem : 

m(z; x) is analytic in z £ C \ F, 



m + {z;x) = m_(z; x) 



e 2i(±z 3 +xz) 



for z£f, 



(3.1) 



_ e -2i(|z 3 +xz) 



m[z\ x) = I + O(^) as z — > 00. 

Here m + (z;x) (resp., m_) is the limit of m(z';x) as z' — > z from the left (resp., right) of the contour F : 
m±(z;x) — lim e |o m ( ;z T The relation between the above Riemann-Hilbert problem and the Painleve II 

equation is the following (see, e.g., 0). If we expand 



. . , mi (a;) nf \ 
m(z; x)=I+ v ' + 0( — 

z y z z 



asz-> 00, 



(3.2) 



G 



we have 

2«(mi(a;))i2 = -2i(mi(x)) 2 i = u{x), (3.3) 
2i(m 1 (x)) 2 2 = -2i(mi(x))n = v(x), (3.4) 



where u(x) and v(x) are dehned in (2.1)-( |2.q ). The above Riemann-Hilbert problem is the special case of 
monodromy data p = —q = 1, r = in the standard family of Painleve II Riemann-Hilbert problems. 
Define 

m M (-iifl;s) to > 0, 

a(af,u>) = < s 3 _ (3.5) 

I — m 2 i(— w; x)e3 tu 2:cu ' w; < 0, 

!mi2(—iw, x) w > 0, 

1 ^ (3.6) 

-mn(-iw; ^ei 1 " 3 - 2 *™ to < 0. 

From the jump condition of the Riemann-Hilbert problem ( |3.l| ), a(x, w) and b(x,w) are continuous at w = 0. 
Indeed since to_u in the upper half plane of C is an analytic continuation of m + in the lower half plane of C, 
a(x,w) and b(x,w) are analytic in w. We have the following properties of a, b. 

Lemma 3.1. For all x,w £ R, we have : 

(i). a{x,w) 1 b{x,w) are real. 

(ii). For each fixed to € R, 



(Hi). 



a(x,w) = I + 0(e- cx3/2 ), x^+oo, (3.7) 

b(x,w) = -ei w3 ~ 2xw (l + 0(e- cx3/2 )), x -> +oo, (3.8) 

a(x,w) ~ -1= e t- 3 -^N 3/2 +N— v^ 2 N 1/2 x^-oo, (3.9) 
v2 

&(*,«,) ~ l e i*>*-#\*\>''+\*\*>-W\*f'> a-^-oo. ( 3 .10) 



lim a(x, iu) = 1, lim b(x,w) = 0, (3-11) 

iu — >-\-oo w — >+oo 



lira a{x. w) = 0, lim 6(x,iu)=0, (3.12) 

— > — oo — ► — oo 

a(x,0) = E{x) 2 , b(x,0) = -E{x) 2 . (3.13) 



(iv). 



9 1 

d~x a{X 


to) 


= u(x)b(x,w), 






(3.14) 


lL b{x 


w) 


— u(x)a(x, w) — 2wb(x, to), 






(3.15) 


d i 


w) 


= 2(u(x)) 2 a(x , w) — \Awu{x) - 


-2u'[x))b(x, 


w), 


(3.16) 




w) 


= (— Awu(x) + 2u(x))a(x, w) - 


f (8u; 2 - 2x 


- 2{u{x)) 2 )b(x,w). 


(3.17) 
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(v). 



a{x,w) = -b{x,~w)e^ w3 - 2xw 1 
b(x,w) = -a(x,-w)ei w3 - 2xw . 



(vi). For each fixed y € R, as w — > — oo, 



a(2?/-\/|H + 4w 2 , w) 
6(2y v / H+4w 2 ,w) 



erf(y), 



f M 3 +4y|H 



3/2 



(3.18) 
(3.19) 



(3.20) 
(3.21) 



Proof. The properties (i)-(v) are consequences of (|3.l|) and Lemma 2.1 in pj. The result (vi) is obtained 



by applying the Deift-Zhou method to the Riemann-Hilbert problem (3.1). The specialty of the scaling x 



2yy\w\ + 4w 2 is related to the fact that the exponent term |z 3 + xz of the anti-diagonal entry in the jump 
matrix of the Riemann-Hilbert problem (3.1) has the critical points at z = ±i-^, which are the stationary 
phase points in the asymptotic analysis (see Recalling z = —iw, we see that x = Aw 2 corresponds to 

one of these stationary phase points. By analyzing m{—iw;x) around this stationary phase point, we obtain 
fii2(— iw, x)e^ w3 ~ 2xw ~ — erf(y) and mn{—iw,x) ~ 1 as w — > — oo with the above x. Similar computation 
appeared in Section 10.3 of [Q, and we omit the detailed computations here. □ 

We now define the following functions with parameters w+, u>_, w. 

Definition 3. For each w + , u>_ € R, when w + + w_ ^ 0, set 

a(x, w+)a{x, W-) — b(x, w+)b(x, W-) 



H(x;w+,w-) = \a{x,w+)a(x,w-)- V ' , , k -i> arU F GU e z . 3.22 

When w + + w_ = 0, we use the l'Hopital's rule (note Lemma [O (v)). Also using Lemma 3.1 (hi), set 



G(x; w) 



lim H(x; w, w-) — a(x, w) Fque(s)- 

_ — ►+oo 



(3.23) 



From Lemma 3.1 (ii), H and G have the limit 1 as x — > oo and has the limit as x — > — oo for each fixed 



w+, W-, w. Also theorem 3.3 below shows that they are limits of distribution functions. Therefore H and G are 



distribution functions. These distribution functions interpolate between the functions Fque, Fqoe and Fo of 
theorem IO] , The following results follow from Lemma 3.1 (hi), (vi). 



Proposition 3.2. For fixed i£R, we have 



H(x; it>_|_, W-) — > < 



Fgue(z), 

FqoeW 2 , 

F , 

0. 



w+, w_ 



-oo 



w + = 0, w_ — > +oo or w_ 
w + = W- = 
w + or u;_ — > — oo, 



-oo 



(3.24) 
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? GUE(a;), w^+oo 
: {Y GO y.(x) 2 . w = (3.25) 

0, w — > — oo. 



j4Zso we /lave /or fixed loii/i u> = — max{— w+, —W-}, as w — > — oo, 



F(2xVkl +4u; 2 ;u; + ,w_) ->• { (3.26) 

erf 2 (x), u; + = w_, 



G(2x^\w\+4w 2 ;w) -> erf(ac). (3.27) 

Now the main theorem in this section is that if we take proper scaling of a± — > 1 in t, we obtain the above 
functions in the limit. 



Theorem 3.3. Set w± by 

We have the followings. 
(i). When < a+ < 1 and vj- € R are fixed, 

L(t) - 2t 



a± = l-^. (3.28) 



lim P — <x =G(a:;u;_). (3.29) 

When vjjf- € R and < a_ < 1 are fixed, the limit is G(x;w+). 
(ii). When w± £ R are fixed, as t — ► oo, 

lim P( ^ ~ 2f < sc ) = ff(x;w+,w_). (3.30) 

A special case of (ii) is when w + — — which corresponds to = 1. In this case, the limiting 

shape of PNG has curvature (see jt3) ) . Hence we obtain a one-parameter family of distribution functions for 
fluctuations of a flat curvature PNG. We have the exact values of the means which include Proposition 2.1 as 
a special case when w = 0. 

Proposition 3.4. We have for each w € R, 

xdH(x; w, -w) = Aw 2 . (3.31) 
Proof. By l'Hopital's rule and Lemma (iv), (v), we have 

H(x; w, -w) = {a(x; w)a(x; —w) - y(x)v(x)} F G tje(z), (3.32) 



9 



where 

y(x) = (2u 2 + x — Aw 2 )a(x; w)a(x; —w) — (u 1 + 2wu)b(x; w)a(x; —w) — (u' — 2wu)a(x; w)b(x; —w). (3.33) 



Then by Lemma 3.1 (iv), (v) again, we obtain 

y'{x) ~ a(x; w)a(x; — w), (3.34) 

which implies that 

H(x;w,-w) = {y'{x) - y(x)v(x)} F GVE (x) = -^-|y(a;) F G UE(a;)|. (3.35) 



Hence as in the proof of Proposition 



lim 

x — >oo 



xH(x;w,—w)— / H(y;w,—w)dy 



lim (xy' - xvy - y) Fque(^) = 4u> 2 , (3.36) 



using the asymptotics of u, v, a, b. □ 

Remark. The function G(x;w) appeared as F M (x;w) in (2.21), (2.22) of [Q. Indeed this function was obtained 
in a different point selection process. Namely in the open square (0, 1) x (0, 1), let S = {(t,t) : < t < 1}, 
the diagonal line, and let 6* — {(t, 1 — t) : < t < 1}, the anti-diagonal line. We select AP(t 2 ) points in 
(0, 1) x (0, l)\(SU5 t ), 2P(at) points on ±), and 2P{(3t) points on <S*\(|, \) such that the resulting point 

configuration is symmetric with respect to both i5 and 8 t . Let L m {t; a, (3) be the length of the longest up/right 
path of a random point configuration in this point selection process. In Q, it is proved that for any fixed j3 > 0, 

L M (t; a, 0) — At < \ ^ fF G u E (x), < a < 1, 

Fgoe^) 2 , ot = 1, 



2tV3 

L M {t;a,(3) -2{a + c l - 1 )t 



^/2{a-a- 1 )t 1 ' 



2 



< x erf(x), a > 1, (3.38) 



^§§^<*) - a.!-* (3.39) 



This result is identical to the limits of L(t) when a_ < 1 is fixed. Comparing this with theorems 1.1 and 3.3 



we see that after scaling, L s (t; a, (3) with any fixed j3 and L(t) with any fixed a_ < 1 have the same statistics 
in the limit t — ► oo. Indeed, one can prove more than that. In finite t, L(2t) with a_ = and L M (t; a+, 0) are 
the same (see Remark at the end of Step 1. in Section ||.) 



4 Lattice directed polymer and exclusion process 

In this section, we analyze a certain lattice directed polymer problem which is closely related to the point 
selection model discussed above. 

Let N = {1,2,3, ■•• } and N* = {0, 1, 2, • • • }. Let < q < 1, a± > be fixed numbers such that a±- s /q < 1. 
We denote by g(q) the geometric distribution with parameter q : 

P(g(q) = k) = (l-q)q k , k = 0, 1, 2, • • • . (4.1) 
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When a random variable w has distribution g(q), we use the notation w ~ g(q)- At each site 6 N* x N*, 
we attach a random variable w(i,j) where 



~ g(q), (t,i)GNxN, (4.2) 

w(i,0) ~ *eN, (4.3) 

w(0,j) ~ fl(a-Vg), jeN, (4.4) 

«;(0,0) = 0. (4.5) 

We call a collection tt of sites in N* x N* an up/right path if when <G 7r, either (i + l,j) € 7r, or + e 7r. 
Let Path(N) be the set of all up/right paths from (0, 0) to (N, N). Define 

X(JV)=max{^] w(i,j):7r£Path(N)}. (4.6) 

The special case a+ = a_ = was introduced by Johansson in P]; the above model adds a special row and 
column to his model. In , it is proved that 

lim J X W-$) N < x ) = Fgue(x) , a+ = a _ = o, (4.7) 



where 

m) = -i f » ct 0? = — i — 7= — • (4-8) 

It is shown in Q that this directed polymer model can be interpreted as a growth model in 2-dimensional 
space, or as a discrete exclusion process. The corresponding exclusion process in our case is the following. We 
use the notation + for the location of a particle. If a site is vacant, we use the notation — . Initially there 
are particles at the sites {• • • ,—4,-3,-2} U {0} C Z. Hence the initial configuration on Z can be written 
as (• • • , +, +, +, — , +,—,—,—)•■■) where the leftmost — is at the site —1 and the rightmost + is at the site 
0. At each (discrete) time, the rightmost particle jumps to its right site with probability 1 — a+y/q 7 and the 
leftmost hole jumps to its left with probability 1 — ct-^/q, while in the 'bulk', each particle jumps to its right 
(equivalently, a hole jumps to its left) with probability 1 — q if its right site is vacant. 

As in the point selection model, when a± are small enough, X(N) would have the limiting distribution 



Fque as in the case when a± = 0. For general a±, we obtain results parallel to those in Theorems 1.1 and 3.3 



After the following changes, we obtain the same limiting results as in Theorems LI and 5J5 
(i). Every limit is taken as N — > oo. 



(ii). The scaled random variable in Theorem 1.1 (i) and Theorem O 

X(N)-^q)N 
cr(g)A^/3 

where fx(q) and o~(q) are defined in 



(4.9) 
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(iii). The scaling of a± in Theorem 3.3 is now 



(iv). The scaling in Theorem 1.1 (ii) is now 



*W Ka<±, (4.11) 



p(q,a)VN V? 
where 

^'"^d-^Ki-^- 1 )' F^FvP) ■ (4 - 12) 

We note that rj(q,a) > n(q) for a > 1. 

Remark 1. In the above model (also in 0]), only the paths ending at the diagonal point (N, N) are considered. 
It is of interest to obtain similar results for path ending at general point (M, N) . The difficulty in this general 
case comes from the fact that we need asymptotics of orthogonal polynomials with respect to the weight function 
(1 + y/qz) M {\ + y/qz- 1 ^ (see Step 2. of Section |). This weight function is not real for \z\ — 1, which makes 
the Riemann-Hilbert method (which we employed to obtain asymptotics) more difficult to analyze. We are 
planning to come back to this problem in the future. 

We note that in |)|. Johansson was able to obtain results for this general case when a+ = a_ = 0. He 
have used different determinant expression (of Fredholm type rather than Toeplitz type) involving different 
orthogonal polynomials (which is discrete), and did not involve the non-real weight function. 

Remark 2. If we take the limit as q — > 1, we obtain exponential random variables instead of geometric distribu- 
tion (see H). In order to investigate the exponential random variables case, we set q = 1 — 4 and I = xL, and 
let L — > co in the determinantal formula ( |5.S| )-( [5TC| ). This double scaling limit is not carried out yet, which we 
plan to do in the future. 

Remark 3. If we take the limit as q — > 0, we have rare events, hence we obtain the Poisson process discussed 
above. Indeed, by setting ^fq = and letting N —> oo, we recover the Theorems [Tl] and |J from (pl|)-([4.11[). 



5 Proofs 



In this section, we sketch how one can obtain Theorems [LlJ and |3.3| using the results of |3j, |J. We split the 
proof into three steps. 

Step 1. We will prove the following formula : 

P(i(t) < = e-^ +a -»- t2 {D[ a+a-D'^) (5.1) 

where 

D[ = EueU(l) det(l + a+U)(l + a^)e 2tRcTl ^ u K (5.2) 
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Suppose a + a_ < 1. In the lattice directed polymer model discussed in Section ^, add a random variable 
w(0,0) ~ g(a+a~) at the site (0,0). Let X + (N) be defined by the formula (4.6) with this new random 
variable added. Now this is a special case of the model □ of (7.5)-(7.7) discussed in Section 7 of || : take 
W = W = {0, 1, • • • ,N} C N*, and take q = a-,<f = a + and q$ = t£ = y/q for 1 < j < N (in |], we have 
taken W, W C N, but simple translation makes no change.) From Theorem 7.1 (7.30) of 0], we obtain 



(5.3) 



P(X+(N) < 1) =(1 - a+a-)(l - a +y /q) N [l - a.^qf (1 - q) N * 

■ ^ueu(i) det{(l + a+E/)(l + a_tft)(l + ^U) N (1 + Jqtf) N }. 

Now we set ^fq — jj, and take the limit N — ► oo. Then inside the square {1, ■ • • , N} 2 , we obtain a Poisson 
process of parameter t 2 : ¥(w = 0) = 1 - -jfa, P(w = 1) = (1 - and P(w >2) = ^. The probability 

that there are k points in the square {1, • • • , iV} 2 is equal to 

N-k, f_ f_\ k ^ e' t2 {t 2 ) k 



N 2 J 



(5.4) 



N 2 'N 2 ) k\ 

disregarding the events of having more than one points at one site which has probability in the limit. Similarly 
we obtain Poisson process with parameter a+t on the bottom edge (not including the origin), and Poisson process 
with parameter a-t on the left edge (not including the origin). At the origin, we have a geometric distribution 
with parameter a + a_. Let L + (t) be the length of the longest up/right path in this process. This is related to 
Lit) by 

L + (t) = L(t)+ X , x~ 5 (a+a-), (5.5) 
since any up/right path will include the g(a+ct-) points in the lower-left corner. From ( |5.3| ), we obtain 



(5.6) 



where D[ is defined in (5.1). 

In order to obtain the formula for L(t), set 



Q(x) = HHt) < l)x l , Q + {x) = HL + (t) < l)t> 



l>0 



l>0 



Then using (5.5), we obtain 

Q+(x) = Y F ( L (t) <l-k)P(x = k) = (l- a+a-)(l - a+a^xy 1 Q(x). 



(5.7) 



(5.8) 



0<fc<2 



Thus by comparing the coefficients, we have ( |5.l| ) for a + a < 1. 

Observe that the first quantity in the right-hand side in (^l]) is entire in a + , a_ , and the second quantity is 



polynomial in a+, a_, hence is entire. Since both sides of (5.1) agree analytically for a + a_ < 1 and the right 



hand side is entire, they agree in general a± where they both converge and are defined. Thus (5.1) holds for 
< a± . Similar consideration yields the following formulae for the probability of the lattice directed polymer 
problem : 



(X(N) </) = (!- <* +y /q)» (1 " a-^r (1 - lY W ~ a+a.T/.J 



(5.9) 
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where 



T! = E ueu(l) det{(l + a+U)(l + a_tft)(l + ^qU) N (1 + y/W*)"}, (5.10) 



and a± are subject to the constraint < a± < l/y^g. 

Remark. It is the formula (^J) that makes connection with the process mentioned in Remark of Section [3| 



In f3j (4.16), it is proved that W(L®(t; a, 0) < 21) is given by ( |5.l| ) with a+ = a and a_ = 0. It is not clear why 
these two processes should be the same. 

Step 2. Let 

D l = E U€U(l) e* t **'*V\ (5.11) 

which is in another form, the / x I Toeplitz determinant det(cj-k)a<j,j<i where Cj is the Fourier coefficient of 
e 2tcos0_ j n |F3 aS y m ptotics of in the double scaling limit as t, I — > oo was studied; hence the second step 



in the proof of Theorem 1.1 is to eliminate the term det(l + a+U)(l + a-W) from the integrand in ( |5.l| ). This 
step is established in Theorem 3.2 of jjj. Let 7r„(2;) = z n + ■ ■■ be the n th monic orthogonal polynomial with 
respect to the weight e t( - z+z ^ ] dz / '{2-kiz) on the unit circle \z\ = 1 where z G C : 

^„(z)^ m (z)e t ^+^ 1 )-^ = (S nm JV n , (5.12) 

|*|=1 

for some constants N n . Note that since the weight function is real, all coefficients of n n are real. Wc define 

7f*(z) = Z n TT(z~ 1 ). (5.13) 

In Theorem 3.2 of ||, using the Weyl integration formula for U(l) and familiar Vandermonde type argument 
together with the relations between orthogonal polynomials on the unit circle and those on the interval, it is 
proved that for a+a_ =/= 1, 

D , = <(-o:+K(-Q=-) - a+a-n(- a +)M- a -) Di ,g ^ 

1 — 

For = 1, l'Hopital's rule applies implying with a + = a, a_ = 1/a, 

£>; = {(1 - JM-aM-eT 1 ) - a^f-aM-a -1 ) - a _1 7rj(-a)7rj(-a _1 )}£)j. (5-15) 
We obtain similar formulae for the lattice directed polymer problem: replace the weig ht e t{z+z 1} by (1 + 

Step 3. The remaining task is to obtain the asymptotics of Di and 7T;(— a) as /, N — ► oo in a proper rate. 
This is obtained in (|, || by applying the Deift-Zhou steepest descent method (see |], ||) to the Riemann-Hilbert 
problem for the orthogonal polynomials tti (z) . These asymptotic results are summarized in Section 5 of [Q . 
Theorem L4 follows by plugging in these asymptotics into ( 5.14 ), ( |5.15 ); we omit the calculations. There are 



similar asymptotic results for the lattice directed polymer problem which yield Theorem 3.3; see Proposition 
3.2 and subsequent remarks in Q. 
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